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ABSTRACT 


The deep ocean is characterized by sound propagation that can support wide-area 
surveillance through the use of distributed acoustic sensors near the seabed. Such a deep- 
water sensor network is potentially enabled by phenomena such as Reliable Acoustic 
Path (RAP) and Deep Sound Channel (DSC) through which undersea network nodes can 
transmit and receive data across long distances. To provide a theoretical understanding of 
system effectiveness, the Bellhop acoustic modeling program is used to predict sensor 


coverage and communications range. 
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I. INTRODUCTION 


Traditional undersea acoustic surveillance systems such as the SOund 
Surveillance System (SOSUS) [1] require expensive and vulnerable electrical cables to 
connect sensors to shore sites where the data can be processed. Surveillance sensors 
could instead use acoustic modems to transport data. An acoustic network can move 
sensor data across long distances through a series of acoustic modems to a gateway node 
such as an Unmanned Surface Vehicle (USV) or surface buoy. The data then can be 
telemetered via satellite to remote locations for further evaluation. Persistent sound 
propagation features such as Reliable Acoustic Path (RAP) and Deep Sound Channel 


(DSC) channel can be exploited in the design of a deep-ocean sensor network. 


This thesis provides a theoretical analysis of the potential for a deep-water 
acoustic network’s effectiveness. The Bellhop acoustic modeling program predicts 
system characteristics such as effective range for various communication frequencies and 


source/receiver depths. 
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Il. BACKGROUND 


Acoustic modems allow operators on the shore to remotely monitor deep ocean 
phenomena. For example, the Deep-ocean Assessment and Reporting of Tsunamis 
(DART) system aids the National Oceanic and Atmospheric Administration’s (NOAA) 
ability to remotely monitor pressure events near the ocean bottom uniquely associated 
with a passing tsunami wave. In addition, the U.S. Navy’s Seaweb system has 
demonstrated communications to and from sensors, Unmanned Undersea Vehicles 
(UUVs), and submarines in shallow water less than 300 meters. Deep-ocean 
environmental characteristics suggest the possibility of extending Seaweb functionality 


beyond the littoral and into the deep sea basins. 
A. DART TSUNAMI WARNING SYSTEM 


Tsunami detection capability gained attention with the December 26, 2004 
tsunami near Indonesia. After this devastating event, which caused the deaths of over 
300,000 people, the enhancement of NOAA capabilities was desired to “forecast the 
impact of tsunamis on coastal areas in time to save lives and protect property” [2] 
worldwide. One of the results was the expansion of the seven-node DART system into 
the 39-node DART II system. The technological advancement and increased coverage 
from DART to DART II help provide “a critical portion of NOAA’s tsunami forecast, 


warning, and mitigation system” [2]. 


The DART II components and operation are illustrated in Figure 1. The system 
begins with the Bottom Pressure Recorders (BPR) scattered on the seabed throughout the 
deep ocean at stations shown in Figure 2. Upon sensing a spike in pressure exceeding 
some detection threshold, the BPR sends a report via acoustic modem to a moored buoy 
on the surface. The buoy transmits the report via Iridium satellite to a Tsunami Warning 
Center (TWC). At the TWC, the detection is further evaluated using the MOST (Method 
of Splitting Tsunamis) method to predict the likelihood and effects of a tsunami 


occurrence. The information is then distributed to the public as a tsunami warning or 


general information. The DART II system was completed in March 2008. 


oad] ude Linda nu tcariie cy aad 





Figure 1. NOAA DART II system [From 3] 





Figure 2. NOAA DART station [From 4] 


The DART II acoustic modems provide for two-way communications between the 
moored buoy and BPR “to allow for two way transmission on demand” for the 
“measurement and reporting of tsunamis with amplitude below the auto-reporting 
threshold” [2]. The amount of data the BPR sends can be controlled to help meet the 
requirement for “very long life with few batteries, and an absolutely guaranteed transfer 
of sensor data into the modem, ready for transmission” [5]. Although the DART II 
“application does not require the transmission of large volumes of data nor does it require 
high data rates” [5], the ability to send only the most crucial information from the bottom 


of the ocean is adequate for tsunami warning. 
B. SHALLOW SEAWEB 


The acoustic modems employed in the DART II system are the product of a U.S. 
Navy Small-Business Innovative Research (SBIR) contract with industry. Using the same 
modems, the Navy has deployed over 50 Seaweb networks around the world in waters up 
to 300 meters deep. The Seaweb system begins with a sensor or undersea vehicle (such as 


a submarine) sending and receiving messages through an acoustic modem. The messages 


) 


are transmitted acoustically over the network of modems to a gateway node at the sea 
surface that can further relay, via satellite, the message to the shore or another ship. Two- 
way communications are available and the modems been shown to “effectively 
communicate with a +30 kt platform operating with a modem designed for a standard 9- 
14 kHz band” [5]. Although Seaweb has the potential for use in the deep ocean, efforts to 


date have been limited to using the system in littoral environments. 





Figure 3. Example Seaweb system [From 6] 


OF DEEP OCEAN CHARACTERISTICS 


The deep ocean is characterized by persistent features conducive to long-range 
acoustic communications. Acoustic modem systems, such as Seaweb, could exploit well- 
known phenomena, such as Reliable Acoustic Path (RAP) and Deep Sound Channel 
(DSC). Acoustic reciprocity in a static environment implies that acoustic propagation is 
bidirectional; thus, most acoustic channel properties apply to both transmitted and 


received signals. 
i. Reliable Acoustic Path (RAP) 


Sound traveling through the deep ocean can be received through various 
propagation paths (such as direct path, bottom bounce, surface ducts, etc.) depending on 


the environment and source/receiver depth. These multipath arrivals are “a common 


occurrence in long-range propagation” [7]. While the accumulation of the different 
arrival paths can cause the received signal to deteriorate, “often one path will be 
dominant, and the transmission loss corresponding to it will be minimum compared to 


other possible paths” [7]. One such dominant path is Reliable Acoustic Path (RAP). 


The RAP phenomenon occurs when a source is placed at a greater depth than a 
receiver in the deep sea. Some possible RAPs are shown in Figure 4. As a RAP is 
“sensitive ... neither to near-surface effects nor to the varying losses on reflection that 
characterize bottom-bounce propagation” [7], it does not suffer scattering losses 
associated with surface interference or the losses caused by bottom reflection and 


absorption. A RAP’s transmission loss is relatively low compared to other paths. 





Figure 4. “Reliable” acoustic paths from a deep source to a shallow receiver in the deep 
sea [From 7] 


The RAP characteristics allow for a bottom-mounted or near-bottom device (such 
as an acoustic modem) to consistently send strong signals to a receiver near the surface. 
Because of acoustic reciprocity, an acoustic sensor at or near the bottom could also detect 
objects near the surface. A scenario where this may be useful is with a bottom-mounted 
sensor detecting a shallow (relative to the bottom of the deep ocean) submarine at 


predictable RAP distances. 


A bottom-mounted sensor uses the RAP to its advantage. A bottom-mounted 
sensor can monitor a large volume of water within the RAP. It is estimated that the 
diameter of the surface coverage is approximately 10 times the water depth. As the noise 
sources outside the RAP volume are strongly attenuated by reflection, absorption, and 


scattering, the dominant signals are from within the detection range. Also, with the low 


ambient noise level at the sea floor, a bottom-mounted device can detect targets or 


receive data within the RAP volume at high SNR (Signal-to-Noise Ratio). 
2. Deep Sound Channel (DSC) 


The Deep Sound Channel (or DSC) is caused by the sound speed gradient shift 
from the surface to the bottom of the ocean. Water temperatures are higher near the 
surface, which causes sound velocity to increase. As depth increases, temperature 
decreases, causing sound velocity to decrease. At some minimal velocity, the sound 
speed gradient goes from negative to positive. At these depths, temperature is isothermal 
and the effects of increasing hydrostatic pressure becomes the dominant influence on 
sound speed, causing it to increase with depth [8]. The depth at which there is minimal 
velocity is called the DSC axis “toward which sound rays are continuously bent by 
refraction” [8]. The change in gradients creates a duct where the sound rays refract back 
and forth between axis upper and lower boundaries (see Figure 5) with low transmission 


loss. Figure 6 shows worldwide DSC axis depths, with a typical axis depth at 1000 m. 
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Figure 5. Deep water sound speed profile with ray trace in DSC, source at 500 m 
[From 9] 








Figure 6. Worldwide DSC axis depths [From 10] 


DSCs have been used to provide long-distance undersea communications. After 
World War II, Ewing and Worzel demonstrated that a small explosive charge could be 
detonated at the DSC axis and sound transmitted over “long ranges ... made possible by 
the natural sound channel which exists in the oceans” [11]. They proposed a system of 
receiver stations that could locate “planes, ships, and life rafts in distress in the open 
oceans” to “within 1 mile” [11]. At the same time, similar DSC experiments were done 
by Brehovskikh [12]. More modern uses of DSCs have been for monitoring nuclear tests 


and for undersea tomography. 


DSCs occur where the sound speed profile has a distinctive deep water shape. The 
Munk canonical equation is “a useful guide to propagation in temperate latitudes” [13]. 
The Munk equation is given in Equations 1 and 2 [After 13] and is plotted in Figure 7 as 


a “typical” deep water sound speed profile. Thus, 


c(z) =1500.0[1.0+ e(e "+7 -1)] (1) 
where 
_ 2(z —1000) Q) 
1000 


and ¢ = 0.00565 and z = water depth. 
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Figure 7. 4000-m Munk sound speed profile 


D. DEEP SEAWEB CONCEPT 


Acoustic modem networks such as Seaweb could exploit low-transmission loss 
paths such as Reliable Acoustic Path (RAP) and Deep Sound Channel (DSC). The 
conceptual system begins with an acoustic sensor at or near the sea floor. A detection of a 
surface contact or submerged contact is made using RAP. The sensor then transmits via 
acoustic modem to a sub-surface repeater node. The node is another acoustic modem that 
can relay the data to an adjacent node using the DSC channel. The data can then either be 
sent to another node or transmitted upward towards a gateway node such as an USV. The 
gateway node then relays the data, via satellite or other radio means, to a shore or ship for 


processing. 
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Deep Seaweb concept 





Figure 8. Proposed Deep Seaweb system [6] 


A Deep Seaweb system has its advantages. Since DSCs are found worldwide, as 
shown in Figure 6, the sub-surface nodes can be placed in a variety of locations providing 
great flexibility to the mission planner. In addition, the use of DSC long-range 
communication enables nodes to not necessarily be close together but positioned for 


optimal area coverage. 


The characteristics of Deep Seaweb make it appealing for military and civilian 
applications. In areas such as shipping lanes, sub-surface nodes on the DSC axis 
(typically 1000 m) are immune to the risk of being hit by a passing ship. For clandestine 
operations, deep-water placement of nodes and use of low-profile USVs minimizes the 
number of objects broaching the surface, reducing the possibility of counter-detection. 
Furthermore, use of the DSC means that the USV does not necessarily need to be near the 
bottom-mounted sensor. The data transfer can occur with long standoff and be 


transmitted to the surface at an area that is appropriate. 
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Deep Seaweb has the potential to be rapidly configurable. In [14], an algorithm 
was developed for discovery of Seaweb nodes following an ad hoc deployment. The 
network-layer routes are initialized according to an optimization cost function. If a node 
is inadvertently removed or damaged, the Seaweb discovery process permits automatic 


healing of the network routes. 
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Hil. ACOUSTIC PROPAGATION MODEL AND LINK MARGIN 


A. BELLHOP 


Bellhop is a Gaussian beam tracing program [15] used in this thesis to model 
transmission loss from an acoustic source. Bellhop was chosen for this analysis as “it has 
proven to be an accurate modeling tool for high-frequency (>1 kHz) transmissions” [15]. 
The Bellhop program was developed by Porter and Bucker at the Space and Naval 
Warfare Systems Center in San Diego in 1987 [15]. 


The AcTUP£ (Acoustic Toolbox User-interface & Post-processor) program, 
created by Maggi and Duncan of Curtin University of Technology in Perth, Australia 
[16], provides the front-end support for the version of Bellhop used in this analysis. The 
principal environmental inputs into AcTUP? include a range-dependent Sound Speed 
Profile (SSP), water depth (z) and water density (p;). The program allows for the 
consistency of the sea floor (such as sound speed cp, attenuation coefficient, and density 
p2) to be added to measure the reflection from the bottom. Acoustic system inputs are 


frequency (f), source depth (zs) and receiver depths (za). 
I Bs Bellhop Validation 


To demonstrate the validity of the Bellhop program at frequencies as high as 10 
kHz, several cases were tested to compare the Bellhop model with the expected 
theoretical results. The MATLAB program in Appendix A calculates the transmission 
loss from a source to a receiver, using an isospeed sound velocity profile. The program 
calculates the amplitude of pressure (P), with respect to depth and range, to find the 
transmission loss (7L). The program uses Lloyd’s mirror (otherwise known as the method 
of images) to find the direct path and both surface and bottom interference pressure 
contributions. These pressures are then combined to calculate total pressure from the 
source to the receiver. The sea surface and sea floor are assumed to be flat and p, is 


assumed to be constant using Kinsler et al.’s value of 1024 kg/m’ [17]. An example of 
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the Lloyd’s mirror geometry is shown at the surface in Figure 9. The geometry at the sea 
floor is similar, except in the opposite direction. Included in Table | are the symbols used 


in the supporting equations below. 


Imaged 
Source sin O 


Sea Surface 
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Source 


Receiver at 
depth zp 


P (rz) 


Sea Floor 


Figure 9. Lloyd’s mirror geometry [After 18] 


Table 1. | Lloyd’s mirror symbols 





Symbol Equation Definitions and Units 





ZR Receiver depth (m) 





Zs Source depth (m) 





Zz Bottom depth (m) 











ig Horizontal range (m) 
R Range from surface ( at r = 0) to receiver (m) 
Ri; Direct path distance from source to receiver (m) 





Ros | Surface reflected distance from source to receiver (m) 
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R23 | Bottom reflected distance from source to receiver (m) 
6 Angle between R and r (rad) 

cA Frequency (kHz) 

a) Frequency (rad/s) 
A Wavelength (m) 

Cc] Sound speed in seawater (m/s) 

C2 Sound speed in sea floor (m/s) 

pl Seawater density (1024 kg/m’) 

p2 Sea floor density (kg/m’) 

vy Intermediate variable used to calculate Re (m*/kg) 
Y2 Intermediate variable used to calculate Re (m*/kg) 
B Intermediate variable used to calculate Re (m'') 
Os Source angle relative to horizontal used to calculate g (rad) 
Q Angle of reflection relative to ground (rad) 
y Complementary angle of g (rad) 

Rr Bottom reflection coefficient (1) 
k Wave number (m”) 
a Attenuation coefficient (dB/km) 
P(r,z) | Pressure as a function of range and depth (wPa) 
P,er | Reference pressure for water (1 Pa) 
TL Transmission loss (dB) 
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To find the pressure at a given range and depth caused by direct path and surface 
interference, first the travel distances R; and R2s are found over range r using Equations 3 


and 4 [After 18]: 
Rar +@,—2;)° (3) 


Rys = r + (Zp Ze). (4) 
Secondly, to account for the attenuation of sound as it travels through seawater, 
the frequency-dependent attenuation coefficient a is estimated, using Equation 5 [After 


18]. Equation 5 is a derivation of the curve in Figure 10. 
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Figure 10. Acoustic attenuation coefficient as a function of frequency [From 18] 
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Next, the wavelength and, thus, wave number are found using Equations 6 and 7. 


ae 


A 6 
f (6) 
20 

aa @ 


The total pressure amplitude is calculated using Equation 8 [After 18]. The 
contribution by surface interference is negative compared to that by the direct path. The 
surface is treated as a pressure release boundary and there is a phase shift, which makes 
the reflection coefficient Rr = -1 [19]. 

22 ape (8) 
R, Rug 

Lastly, the transmission loss 7L is calculated by Equation 9. 

TL = 20 log(—) (9) 
ref 

For the bottom 7L calculation, the equations are the similar except that R2g 
replaces R2sto signify the change in the reference plane for Lloyd’s mirror. Also since the 
sea floor is a fluid-fluid boundary, the reflection coefficient Rr must be found to 
determine how much of the signal is reflected from the sea floor toward the receiver. The 
reflection coefficient is found by comparing the density of seawater to that of the bottom, 
as well as the sound speed change [After 18]. 

R, = PHP (10) 
Pal + PAY 


The variables y; and y2 are defined as: 


(2) 
y, =—cos(y) (11) 

C, 

Cc 2 
= u 12 
‘ cy sin’ (y)-1 ow 
where 

B=" siny) (13) 

Cc 


and 


1 
=—-@. 14 
ae (14) 
Furthermore, the angle g is defined as: 
p= tan“!(25* 7) = tan“!(25* 72) (15) 
Ltr r 


Sea Surface 


r 


Receiver at 
depth zp 


P (rz) 





Sea Floor 


Figure 11. Definition of 0s, g, and y used to calculate bottom reflection 
coefficient Rp 


With Rr known, Equation 8 is calculated as before and the transmission loss 7L is 


calculated by using Equation 9. 
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Two examples are used to compare the TL calculated by the above equations with 
that of the Bellhop model. The first example compares the TL caused by surface and 
direct path contributions. A 10-kHz signal emanates from a source zs = 20 m and is 
received at ze = 200 m. The environment consists of an isospeed sound profile with c; = 
1500 m/s, p; = 1024 kg/m’, and water depth z = 5000 m. The two results in Figure 12 are 


similar. 
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Figure 12. Comparison between direct path and surface 7Z contributions using 
theoretical and Bellhop values 
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The second example shows the bottom reflection and direct path contributions. A 
10-kHz signal is now projected from a source zs = 4975m and received at zr = 4800 m. 
The water column has an isospeed sound profile where c; = 1500 m/s. The sea floor is at 
a depth of z = 5000 m has properties of c2 = 1800 m/s and p>: = 1843 kg/m*. The 
theoretical and Bellhop 7Z closely track each other in Figure 13. 
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Figure 13. Comparison between direct path and bottom 7Z contributions using 
theoretical and Bellhop values 


The Bellhop program 7Z values match the theoretical values given the same 


initial conditions. Bellhop is considered validated for the remainder of this study. 
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B. SONAR EQUATION 


The ability of a receiver to detect a signal requires it to be above the 
environmental and equipment noise that can limit detection. The passive sonar equation 
in Equation 16 gauges the ability of a transducer to detect acoustic signals. Rearranging 
the passive sonar equation [After 7], the required level needed to detect a sound is given 
by: 

SNR = SE— NL+ DI -TL (16) 


Table 2.__ Passive sonar equation symbols 





Symbol Equation Definitions and Units 





SNR_ | Signal to Noise Ratio (dB) 





SL Source level (dB re 1uPa @ | m) 





NL Noise level (dB) 


DI Directivity index (dB) 





TL Transmission loss (dB) 














Each component of the passive sonar equation contributes positively or negatively 
to the Signal-to-Noise Ratio (SNR). The components are measured in dB, requiring them 
to be normalized by a reference value. A_ typical dB calculation 


Value 


is 20 log( : 
Reference Value 





The source level (SZ) is the radiated pressure level relative to 1 wPa measured at 1 
m away from the source. The SZ is the standard definition for the amount of energy the 
source puts in the water. Noise Level (NZ) is the additional received energy caused by 
outside environmental events (such as wind, shipping, or biological). A higher NZ makes 
the signal harder to detect and detracts from SNR. Traditionally NZ is estimated using 


frequency-dependent Wenz curves which account for sea state, shipping, and wind. 
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Figure 14. Wenz curves to determine ambient noise level [From 20] 


The directivity index (DJ) is the ability of the transducer to physically or 
electronically direct its received response in space. An omni-directional transducer can 
detect a signal in all directions and has DJ = 0. A directional transducer detects better in 


a certain bearing and has DJ > 0 dB. 


Lastly the transmission loss (7L) is the loss suffered while the signal travels 
through the medium. The 7Z depends on water column, water depth, and sea floor 
characteristics. Again because of reciprocity, two identical transducers being used as a 


source and receiver are interchangeable in terms of 77. 


For a given sonar system in a given noise environment, as long as the SNR 


exceeds the detection threshold, the signal can be detected. The use of an acoustic model 
22 


such as Bellhop can identify sensitivities caused by changes in environmental factors. 
Using the same initial conditions as for Figure 13, a hypothetical assumption is made that 
the value of SL—-NL+DI is 50 dB. For a detection to occur, TZ needs to be < 50 dB. For 
TL values higher than 50 dB in Figure 15, the losses are too great and detection will not 


occur. 
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Figure 15. Detection range limit as a function of 7Z for surface interference example 


23 


THIS PAGE INTENTIONALLY LEFT BLANK 


24 


IV. RELIABLE ACOUSTIC PATH (RAP) 


Reliable Acoustic Path (RAP) in the deep ocean can be exploited for sensing and 
communications. RAP characteristics vary with factors such as transducer placement, 
water depth and signal frequency. A parametric analysis can reveal how sonar 
performance in a RAP environment is affected by these factors. The following studies use 
the property of reciprocity, that is the TL experienced from a source to a receiver is the 


same as if the transducer roles were reversed. 


An acoustic modem and sensor apparatus is anchored at the bottom of the ocean, 
as shown in Figure 16. It is assumed that the acoustic sensor is a transducer baffled such 
that it provides a hemispherical response. The hemispherical response limits the 
transducer’s ability to detect those signals arriving from 0° to 90° above the horizontal 


axis. 


Acoustic modem 







transducer 
5m 
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Float 
Electronics 
Sm 


Anchor with 
batteries 


Sea floor 


Figure 16. Undersea modem and sensor apparatus for RAP analysis 
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A. RECIPROCITY 


For a 4000-m water column, a Munk profile with the same values as Equation | 
and similar to Figure 7 is considered. The bottom is assumed to be gravel with sound 
speed cz = 1800 m/s and bottom compression attenuation coefficient a) = 0.6 dB/A. Shear 
losses in the water column and bottom are neglected. A 100-Hz source radiates from a 


near-surface depth of 30 m, representing a surveillance contact. 


Transmission loss, as a function of range and depth, is shown in Figure 16. The 
brighter areas represent lower transmission loss and the RAP “region” is visible. At 


approximately 30 km, the RAP ceases to converge with the bottom. 
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Figure 17. TL of a 100-Hz, 30-m source using a 4000-m Munk sound speed profile 


The environment is assumed to be static and acoustic reciprocity applies. To 
demonstrate the property of acoustic reciprocity, a 100-Hz source is placed at a near- 
bottom depth of 3995 m. Figure 18 shows that the RAP reaches the surface to ranges of 
30 km, the same range as in Figure 17. The RAP region seen in Figure 18 exhibits the 


radial cross-section of the classic “tea cup” surveillance volume for a bottom sensor. 
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Figure 18. TL of a 100-Hz, 3995-m source using a 4000-m Munk sound speed profile 


B. DEPENDENCE ON OCEAN DEPTH 


With a sensor height of 10 m above the sea floor, the ocean depth is now varied 
from 500 m to 4000 m and Bellhop 7Z plots are produced for selected ocean depths. The 
range of maximum signal level (Rr4p) is plotted versus ocean depth in Figure 19. The 
Rrap increases as bottom depth z increases. The plot has three distinct regimes where 
linear relationships can determine Rayp as a function of bottom depth. From a depth of 0— 
3.25 km, Rrap ~ 5.8 * z. Deeper still from 3.25—3.5 km , Rraap ~ 6.7 * z. From 3.5—4 km, 
Rrap ~ 7.6 * z. As the bottom depth gets deeper, the gradient for Rayp gets larger. The 
diameter of the surveillance area is 12-14 times the bottom depth which validates the rule 
of thumb stated in Chapter II. For example, a single sensor deployed in 4 km water will 


have a surveillance area of 2800 km?! 


ai 














Figure 19. Ocean depth vs. Reap, 100-Hz, 3995-m source using a 4000-m Munk 
sound speed profile 


C. SENSITIVITY TO MIXED LAYER VARIATIONS 


The upper portion of the water column is subject to physical mixing by wind and 
surface waves. The sound speed in the mixed layer tends to be uniform, with little or no 
vertical gradient. Using a representative temperate latitude SSP (15.5°N, 179.5°E) from 
[13], a mixed layer with depths of 50 m and 100 m is introduced in Figures 20 and 21. 
The remainder of the sound speed profile remains as a 4000-m Munk profile. A mixed 
layer near the surface of the water column does not appreciably change the detection 


ranges on the bottom of the ocean, as shown in Figures 22 and 23. 
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Figure 20. Temperate latitude 4000-m Munk sound speed profile with 50-m mixed 
layer depth 





Figure 21. Temperate latitude 4000-m Munk sound speed profile with 100-m mixed 
layer depth 
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Figure 22. Temperate latitude TL of a 100-Hz, 3995-m source in a 4000-m Munk 
sound speed profile with 50-m mixed layer depth 
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Figure 23. Temperate latitude TL of a 100-Hz, 30-m source in a 3995-m Munk sound 
speed profile with 100-m mixed layer depth 
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D. FREQUENCY DEPENDENCE 


Using a bottom depth of 4000 m and a source depth of 30 m, the frequency is now 
varied. Figures 24 and 25 show the 7Z for 1 kHz and 10 kHz, respectively. As the 
frequency increases, so does TL due to the frequency-dependent attenuation coefficient a 
from Equation 5. At 10 kHz, the 7L is so great that the RAP does not reach the bottom. 
The significance of this is that acoustic modems that operate in the 9-14 kHz band will 


have to be placed almost directly beneath a gateway buoy to transmit to the surface. 
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Figure 24. TL of a 1-kHz, 3995-m source in a 4000-m Munk sound speed profile 
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Figure 25. TL of a 10-kHz, 3995-m source in a 4000-m Munk sound speed profile 


E. LINK BUDGET 


Another tool that can be used to gauge the effectiveness of signal transmission 
through a medium is link budget analysis. A link budget analysis quantifies the ability of 
the sensor system to detect a signal against a background of noise. Using the passive 
sonar equation from Equation 16, the SNR can be estimated for a sensor’s ability to detect 
the 100-Hz target at depth of 30 m in sea state 3 with normal shipping. From the Wenz 
curves in Figure 14, the NZ is estimated to be 43 dB. The SZ is assumed to be 130 dB (re 
1 uPa @ | m) and the transducer DJ is 3 dB. Equation 16 becomes: 


SNR = 130-43+3-TL (17) 
SNR = 90-TL (18) 
Detection is likely as long as TL is less than 90 dB. Figure 26 shows the detection 


ranges for the above scenario for a given detection threshold. 
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Figure 26. Transmission range of a 100-Hz acoustic modem, 3995-m source and 30- 
receiver in a 4000-m Munk sound speed profile 


F. IMPULSE RESPONSE 


The impulse response was extracted from Bellhop using methods discussed in 
[15]. The impulse response is an indication of the channel time spread produced by 
multipath propagation from source to receiver. The channel time spread is an important 
consideration for communications inter-symbol interference (ISI). Some sample paths are 


shown in Figure 27. Representative impulse responses are shown in Figures 28 and 29. 
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Figure 27. Sample multipath propagation 
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Figure 28. Impulse response for a 7-kHz acoustic modem, 3990-m source and 30-m 
receiver with launch angles of +89° in a 4000-m Munk sound speed profile with 
source-to-receiver range of 500 m 





Time (s) 


Figure 29. Impulse response for a 7-kHz acoustic modem, 3990-m source and 1000- 
m receiver with launch angles of +89° in a 4000-m Munk sound speed profile with 
source-to-receiver of 500 m 
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V. DEEP SOUND CHANNEL (DSC) 


The system considered in the DSC analysis is an acoustic modem tethered to an 
anchor on the sea floor as depicted in Figure 30. The water column is characterized by a 
4000-m Munk SSP (Figure 7), and the bottom characteristics are the same as for the RAP 
analysis. The length of the tether is 3000 m optimally placing the acoustic modem 
placement directly on the sound channel axis at a depth of 1000 m. The standard Seaweb 
modem operates in the 9-14 kHz band and it is initially assumed for this analysis that the 
operating frequency is 10 kHz. 
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Figure 30. Undersea modem apparatus used for DSC analysis 


A. VERTICAL LAUNCH ANGLES 


The acoustic modem vertical launch angles dictate the initial trajectory of 
propagation into the water medium. If the vertical beam pattern is too wide, energy will 
leak out of the DSC, reach the surface (or bottom) and be wasted. A similar sound speed 


profile with the source placed at a 1000-m DSC axis is shown in Figure 31. 
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Figure 31. Example of launch angles from source located in DSC axis at a temperate 
latitude [From 12] 


To find the launch angles at which the rays will stay within the DSC, a method 
developed by Munk [20] is used. The function @ is defined as: 


G=FY(C-C)- eC, (19) 





where 
Cc 
C= 20 
' cos0 m 
and C is the depth at the DSC axis and @ is the launch angle. 
Equation 20 can be rewritten as 
1 
@ =cos '(——) (21) 


l+¢é 
Munk defines the sound channel asg¢=+1. With ¢ =0.00565 as before, 6 = +6°. 


Figure 32 shows the 7L for a 10-kHz source at 1000 m with launch angles of +6°. The 


maximum range is about 25 km on the axis, where 7Z is 100 dB. 
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Figure 32. TL for a 10-kHz, 1000-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile 


B. SENSITIVITY TO PLACEMENT AWAY FROM THE DSC AXIS 


Outside forces, such as currents, can cause the acoustic modem to deviate from its 
intended depth on the DSC axis as shown in [22]. In Figures 33 and 34, the modem 
depth is below the axis by 200 m and 400 m, respectively. When the modem is moved off 
axis, and maintaining the launch angles at +6°, the transmitted signal departs the DSC 
boundaries and the effective range is degraded. At 1200 m and 1400 m depth, the 


maximum range is about 24 km. 
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Figure 33. TL for a 10-kHz, 1200-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile. 
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Figure 34. TL for a 10-kHz, 1400-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile. 
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A similar analysis with the modem placed at depths of 600 m and 800 m shows 
the range decreasing to about 21 km. In addition, Figures 35 and 36 predict range regimes 
wherein there exist shadow zones at the DSC axis when the source is above axis. Placing 
the source directly on the DSC axis is the optimum location for long-range transmission. 
Negligible degradation occurs with the source below the axis. Substantial degradation 


occurs with the source above the axis. 
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Figure 35. TL for a 10-kHz, 600-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile. 
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Figure 36. TL for a 10-kHz, 800-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile. 


C. FREQUENCY DEPENDENCE 


The acoustic modem is repositioned on the DSC axis for optimum performance. 
The frequency is now varied to study the effect on 7L. Figures 37-42 show a progression 


as frequency increases. The increasing attenuation results in less achievable range. 
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Figure 37. TL for a 3-kHz, 1000-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile 





Figure 38. TL for a 5-kHz, 1000-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile 
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Figure 39. TL for a 7-kHz, 1000-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile 
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Figure 40. TL for a 9-kHz, 1000-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile 
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Figure 41. TL for a 14-kHz, 1000-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile 
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Figure 42. TL for a 20-kHz, 1000-m source with launch angles of +6° in a 4000-m 
Munk sound speed profile 
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BELLI 


For select frequencies between 5—13 kHz, the 7Z is plotted versus range. Figure 
43 shows that, as frequency increases, the TL level increases 


HOP. Muri 
T 
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Figure 43. 





TL vs. frequency for a 1000-m source and 1000-receiver with launch 
angles of +6° in a 4000-m Munk sound speed profile 
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D. LINK BUDGET 


The link budget analysis using Equation 18 is done for a 7-kHz acoustic modem 
at a depth of 1000 m. The SZ is 130 dB with a sea-state 3 NL of 18 dB. The DI of the 


modem is 6 dB. The acoustic modem transmits to another modem on the axis. 








SNR = 130-18+6-T7L (22) 
SNR =118-TL (23) 
’ l ly 
$F oe fe 
i 
Figure 44. Transmission range of a 7-kHz acoustic modem, 1000-m source and 1000- 


receiver with launch angles of +6° in a 4000-m Munk sound speed profile 
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E. IMPULSE RESPONSE 


The modeled impulse response for a 7-kHz acoustic modem at the DSC axis is 


given in Figure 45. 


Figure 45. Impulse response for a 7-kHz acoustic modem, 1000-m source and 1000- 
m receiver with launch angles of +10° in a 4000-m Munk sound speed profile with 
source-to-receiver range of 20 km 
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VI. ACOUSTIC UPLINK TO SURFACE GATEWAY NODE 


The final link of the Deep Seaweb network is the communications from the 
subsea network to a surface gateway node, such as a USV or moored buoy. The link 
could come either from the bottom-mounted sensor or from a network node placed in the 
DSC. Figures 24 and 25 show the frequency variations of 7Z from a deep source to a 
shallow receiver. A higher-frequency (10 kHz) modem at the DSC axis would require the 
gateway node to be nearly overhead, a lower-frequency modem allows for RAP usage 
with much greater area coverage. Figures 46-49 show that, as the frequency is reduced, 
the RAP begins to reveal itself. At 7 kHz, the RAP seems fully formed. In conjunction 
with lowering the frequency, increasing the upward DI of the modem can mitigate TL 


effects. 
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Figure 46. TL for a 9-kHz, 3990-m source in a 4000-m Munk sound speed profile 
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Figure 48. TL for a 5-kHz, 3990-m source in a 4000-m Munk sound speed profile 
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Figure 49. TL for a 3-kHz, 3990-m source in a 4000-m Munk sound speed profile 


When the source is a modem placed at the DSC axis, the signal has less distance 
to travel and, thus, suffers less attenuation. Figures 50-53 show that as frequency 
increases, so does TL. However, with less distance traveled, there is less chance for 
spreading before a boundary interaction and less effective range. Where the effective 
range for the bottom-mounted modem is around 25 km, the range for a DSC axis node is 


around 7 km. 
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Figure 50. TL for a 9-kHz, 1000-m source in a 4000-m Munk sound speed profile 
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Figure 51. TL for a 7-kHz, 1000-m source in a 4000-m Munk sound speed profile 
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Figure 52. TL for a 5-kHz, 1000-m source in a 4000-m Munk sound speed profile 
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Figure 53. TL for a 3-kHz, 1000-m source in a 4000-m Munk sound speed profile 
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VI. CONCLUSION 


Existing networks, such as DART and Seaweb, have been used for wireless 
transmission of data underwater. Both have been able to relay data from a bottom- 
mounted sensor over long distances to a surface buoy. Both systems subsequently 
transmit the data via satellite to a shore-based site for further evaluation. A system such 
as Deep Seaweb combines attributes of these systems, and exploits RAP and DSC for 


military and civil long-range acoustic communications. 


The Bellhop acoustic model provides a theoretical understanding of a deep-water 
acoustic network’s effectiveness with variations of depth, sound-speed profile, and 
frequency. Given predictions of effective acoustic ranges, node placement can be 
optimized so an efficient number of nodes can be used for area surveillance and 
communications. The optimal placement of a sensor trying to exploit RAP is on or near 
the sea floor. In the DSC, the optimal node placement is at or just below the axis. The 
gateway node can connect either to a DSC node or to a seabed node. The operating 
frequencies of the acoustic modems should be as low as practical, while still retaining 
adequate spectral bandwidth. Recommend Deep Seaweb be developed with a 


communications band below 8 kHz, e.g., 5—7.5 kHz. 
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APPENDIX A. LLOYD’S MIRROR MATLAB CODE 
clear all 


r =[1:5000]; “receiver ranges (m) 

f= 10; “frequency (kHz) 

ZS = 25; Y%source depth (m) 

zr = 200; “receiver depth (m) 

zbottom = 5000; “obottom depth (m) 

cl = 1500; %isovelocity sound speed (m/s) 
c2 = 1800; %bottom sound speed 

rhol = 1024;% density of seawater (kg/m‘3) 
rho2 = 1843;% density of bottom (kg/m‘3) 
fl = f*10°3; “frequency (Hz) 


w = 2*pi*fl; Yofrequency (rad/s) 

lambda = cl/fl; % wavelength in seawater(m) 

k = 2*pi/lambda; Yowave number 

pref = 1; Yoreference pressure (uPa) 

zsb = zbottom - zs; “ochanges reference point for bottom direct path 
zrb = zbottom - zr; Yochanges reference point for bottom reflection 


thetal = atan((zsbt+zrb)./r); 
theta = pi/2 - thetal; 
thetadeg = thetal * 180/pi; 


gammal = (w/ cl) * cos(theta); 
B = (w/cl) * sin(theta); 
gamma2 = B .* sqrt(cl’2 ./ (c2“2 * (sin(theta).%2)) - 1); 


R = (tho2 * gammal - rhol * gammaz2) ./ (rho2 * gammal + rhol * gamma2); “%bottom 
reflection coefficient 
Rphase = angle(R) ./ pi; 


alpha = 3.3*104-3 + [(0.11*f2)/(1+f%2)] + [(44*f%2)/(4100+f*2)] + £42*3.0*10%-4; 
‘%attenuation coefficient (db/km) 

alphal = alpha/1000; “cattenuation coefficient (db/m) 

alpha2 = alphal/8.7; Yoattenuation coefficient conversion (Np/m) 


Rls = sqrt(r.42 + (zr-zs)*2); % surface direct path distance (m) 
R2s = sqrt(r.42 + (zrt+zs)*2); % surface reflected path distance (m) 
pls = [exp(11*k*R1s).*exp(-alpha2*R1s)]./R1s; “pressure caused by surface direct path 
(uPa) 
p2s = [exp(11*k*R2s).*exp(-alpha2*R2s)]./R2s; “opressure caused by surface reflected 
path (uPa) 
ps = pls-p2s; Yototal pressure (uPa) 
a7 


TLs = -20*log10(ps/pref); Yotransmission loss with surface interaction (db) 


Rlb = sqrt(r.*2 + (zrb-zsb)*2); % bottom direct path distance (m) 

R2b = sqrt(r.*2 + (zrb+zsb)*2); % bottom reflected path distance (m) 

plb = [exp(1i*k*R1b).*exp(-alpha2*R1b)]./R1b; “pressure caused by bottom direct path 
(uPa) 

p2b = [R.*exp(1i*k*R2b).*exp(-alpha2*R2b)]./R2b; Yopressure caused by bottom 
reflected path (uPa) 

pb = plb+p2b; “ototal pressure (uPa) 

TLb = -20*log10(pb/pref); “otransmission loss with bottom interaction (db) 


“Yonote pls and plb should be the same since they are both the direct path 
“pressures from the source to the receiver 


figure(1) 

plot (r,TLs,'r'); 

title (‘Surface Interference’) 
xlabel (‘Range (m)'); 

ylabel (‘Transmission Loss (dB)'); 
legend Bellhop Theoretical 


figure(2) 

plot (r,TLb,'k’); 

title (‘Bottom Interfernce’) 

xlabel (‘Range (m)'); 

ylabel (‘Transmission Loss (dB)'); 
legend Bellhop Theoretical 
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APPENDIX B. MODIFIED IMPULSE RESPONSE MATLAB CODE 


fid=fopen('RAPIR.arr'); 

for 1=1:6 
tline=fgetl(fid); 

end 

xy=textscan(fid,' “of %of Y%f Yof Yof Yod %od'); 

x=xy {3}; 

yay 

stem(x,y,'k.') 

xlabel('Time (s)') 

ylabel(‘Amplitude') 

title(Impulse Response’) 
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